using Printf, ProgressMeter
res5 = post1(20)
plt.plot(res5[:,1], res5[:,2], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,3], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,4], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,5], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,6], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,7], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,3] .+ res5[:,6], color=:black, linewidth=1.5, linestyle="-", marker="o", markersize=0);
plt.show()
plt.plot(res5[:,1], res5[:,4] + res5[:,7], color=:black, linewidth=1.0, linestyle="-", marker="o", markersize=0);
plt.show()
np.savetxt("IntForceData.txt", res5)
#=
fout = open("data ndp=6.txt","w")

@showprogress 0.1 "Computing... " for i in [2, 4, 8, 16, 32, 64, 128]
    global nels = i
    include("main1.jl")
    @printf(fout, "%5d\t%.6e\t%.6e\t%.6e\t%.6e\t%.6e\t%.6e\n",
    neq, res5[1,2],res5[1,3],res5[1,4],res5[1,5],res5[1,6],res5[1,7])
end

close(fout)
=#
